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Spin-magnetophonon level splitting in a quantum well made of a semimagnetic wide gap semi- 
conductor is considered. The semimagnetic semiconductors are characterized by a large effective g 
factor. The resonance conditions hui^o = j-iBgB for the spin flip between two Zeeman levels due 
to interaction with longitudinal optical phonons can be achieved sweeping magnetic field B. This 
condition is studied in quantum wells. It is shown that it leads to a level splitting that is dependent 
on the electron-phonon coupling strength as well as on the spin-orbit interaction in this structure. 

We treat in detail the Rashba model for the spin-orbit interaction assuming that the quantum well 
lacks inversion symmetry and briefly discuss other models. The resonant transmission and reflection 
of light by the well is suggested as a suitable experimental probe of the level splitting. 
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The resonance coupling of Landau levels with longitudinal optical phonons (magnetophonon resonance) was the- 
~5 ■ oretically predicted in |l| in the magnetoresistance investigation. The resonance takes place every time when the 
optical phonon frequency is the cyclotron frequency of an electron times some small integer. Thus a possibility has 
been pointed out for an internal resonance in solids. This phenomenon since its prediction has been observed in many 
experiments — see for instance the review Q- 

A possibility of spin-flip transitions of electrons interacting with optical phonons between the Landau levels of 
opposite spin orientations that may be called spin-magnetophonon resonance (SMPR) — has been indicated and 
discussed in a number of papers (see [3, EL la @])- The purpose of the present paper is to discuss the peculiarities 
of SMPR in semimagnetic semiconductors where due to large effective g-factors the corresponding interlevel spacing 
Q ' may be particularly large and therefore SMPR is well pronounced. The condition for the spin resonance has the form 

gii B B = Hujlo- (1) 
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t-H ■ Here fiB is the Bohr magneton, g is the carrier effective g-factor while B is the external magnetic field. 

Many remarkable magnetooptical properties of wide gap semimagnetic semiconductors such as giant exchange 

splitting of the free exciton [7| , giant Faraday effect 0, [& Q , etc are determined by a large splitting of conduction and 

\j2 ' valence bands in magnetic field. This is a consequence of the exchange interaction of band carriers with the electrons 

—i. of the half filled d shell of the Mn ions. In the present paper we will treat as an example the compound Cdi-^Mn^Te 

where the width of the gap between the top of valence band and the bottom of conduction band in the absence of 

magnetic field is given by E g = 1.595 + 1.592a; eV. 

In the presence of an external magnetic field the Mn ion spins are aligned along the magnetic field. Through the 
f^ \ exchange of the Heisenberg type these spins interact with spins of the band carriers. Eventually, in the mean field 
!• ■ model the band carrier dynamics can be described incorporating the exchange interaction only into the enhanced 
. £^ | g-factor. 

There are two competing mechanisms determining the sign and value of the exchange constant (and of the g- 

S-h ' factor) [lfj, llll Il2| . The first mechanism originating from direct exchange interaction between the band and d 

electrons is relatively weak and ferromagnetic. The second one is due to hybridization of d orbitals and band states. 

The latter turns out to be antiferromagnetic and is negligible for the conduction band while for the valence band it 

determines the exchange constant. 

The resonance coupling of Landau levels with optical phonons manifests itself also in a different way though the 
underlying physics is the same. It leads to magneto-optical anomalies both in bulk [13J and in two-dimensionally 
confined systems [TJL [Ij|. Primary concern of Refs. [lj, [l5| was investigation of magneto-optical anomalies of opti- 
cal phenomena in conventional GaAs based heterostructures. It was shown that magneto-optical anomalies in two 
dimensions provide a powerful tool for the electron-phonon coupling investigation in these structures. It was found 
that under the resonance condition with respect to electron-phonon interaction the relevant cyclotron peak splits into 
a doublet. This effect leads to anomalies in optical absorption and reflection (as well as in other optical effects such 
as for instance Raman scattering). 

In the present paper we investigate this effect associated with SMPR, i. e. the magnetophonon resonances due to 
the spin flips. These electron-phonon resonance conditions can occur both for the valence and conduction electrons. 
The exchange constants for the conduction and valence electrons turn out to be different [16| • Though the resonance 
condition leading to the level splitting occurs first in the valence band we will show that the splitting itself is smaller 
for the valence band states than for the states in the conduction band. 

In the next section we consider the level splitting as a formal quantum mechanical problem. This phenomenon 
can be understood in terms of degeneracy lifting of two degenerate states. The energy degeneracy of an electron 



state 2 and an electron in a state 1 plus an optical phonon (see FigflJ is lifted by the electron-phonon interaction. 
We will obtain an expression for the level splitting not specifying the states involved in the relevant transitions. 
In Sec. Mil we determine the states and the energy levels of the conduction electrons taking into account the spin- 
orbit interaction in the Rashba model. This allows to express the level splitting explicitly. We give the required 
estimations in the end of this section. As a possible experimental probe of this splitting phenomenon we propose 
the resonant reflection (transmission) of the light by a quantum well in the Faraday configuration. We consider the 
wave reflection (transmission) due to direct interband transitions, therefore in Sec. II VI we give explicit expressions for 
the wave functions and energies of the valence band states. In section [V] these wave functions are used to determine 
the reflection and transmission coefficient of the light exciting interband transitions in the quantum well. In Sec. IVII 
we discuss applicability of the perturabation theory for solution of our problem. We present conclusive remarks in 
Sec. ETTJ 

II. LEVEL SPLITTING 

We begin with treatment of a formal problem: we will consider two states 1 and 2 and find the self energy of an 
electron in the state 2 due to interaction of the electron with optical phonons. Suppose the energy of the state under 
consideration e 2 is close to E\ + fii^ho (i- e - the electron state 2 and the electron state 1 plus the optical phonon with 
frequency wlo & r e degenerate). This allows us to put aside all other possible electron states. 

Generally, a single quantum well brings about new phonon (vibrational) modes. There could be three types of 
phonons associated with a quantum well [17] : phonons not penetrating into the quantum well, phonons peaking at 
the interface and decaying both in the well and in the barriers (interface phonons) , and phonons confined to the well. 
The phonon Green function in the Matsubara technique can be written as 



v- ^ ,2 / « - ~ri a (zW) e-^ T ^{z)r) a {z') 
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D(v u z,z'^ k ) = -^\C a \ 2 [ ■ "IY - , "y'°v~' , (2) 



where rj a (z) describes spatial distribution of the phonon a branch in the direction perpendicular to the well plane [z 
axis), 0Jk = 2irkT (k = 0, ± 1 . . .) are the Matsubara boson frequencies, |C Q | 2 is the electron-phonon coupling strength. 
T is the temperature; we will use for it the energy units setting fee = 1. 

The electron-phonon interaction with longitudinal optical phonons can be treated in the bulk Frohlich [l8|] model. 
According to the model r/ a (z) — > e lq * z , \C a \ 2 —* 2-Ke 2 hu>i J o/q 2 e* . This approximation in a relatively wide wells 
can be justified noting that the interaction with the interface phonons in this case can be neglected, interaction 
with the confined phonons qualitatively leads to the same result. Therefore, further on we will work in the Frohlich 
approximation. We consider the dispersionless optical phonons, cjlo being their frequency and 

i-i-i. (3) 

where eoo(eo) is the high frequency (static) limit of the dielectric susceptibility. 

The electron self energy in the first approximation of the perturbation theory with respect to the electron-phonon 
interaction can be written as (see the diagram (a) in Fig.©) 
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S2(«£ n ) = —T > — r 5 7- nr (4) 
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where 



and 
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e n = 7r(2n + l)T. (6) 



To calculate the sum over k in Eq.([4]) we use 



and 



G{ie n ) = - / de ^- (7) 
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We define the phonon Green function D R as 



UJ — hidLo + iO UJ + hidLo + iO 






The sum can be represented as 

k J -°° 

Here we have used the identities 

T"£^—=n F (x)-l, T^2^—=-[n B (y) + l], (11) 

*-f ie s - x ^ tLu k -y 

where n F (%) \p-B {y)] is Fermi [Bose] function. The final result is 

2iruj h oe 2 h f n F (ei) - n B {uLo) - 1 n F (si) + n B (uj L o) \ 



^ 2 {ie n ) = F 21 <- - r — — — \. (12) 

e L ie n - (ei — fjt) — «^lo «£« - (ei - Mi + "- w lo J 

Restricting ourselves with the low temperature case T <C ft^LO an d assuming that the state E\ is empty we get 

A 2 /4 
ie n - (ei - n) - 7kj LO 



where 



A 2 = 8 "" Loe2 ^ 21 . (14) 
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For the electron Green function we get with this self energy 

G2{l£n) = ie n -e 2 + v- WWJWn -ex- fiw LO + V) ' (15) 

We perform analytical continuation replacing i£„ — ► e + iO and get for the retarded Green function 

G«( £ ) = 7 e- £l -hu; L o+» 

2W (£-£++^ + i0)(£-£_ +^ + i0) v ; 



where 



£2 + £l + S-WLO 



e± = 7 ± V((£2-£i-^ LO )/2) 2 + (A/2)2. (17) 



As is seen from Eq. ([T6|) we have gotten two poles of the Green function; the level £2 is split into a doublet with the 
energies e±, the spacing between the poles being equal A. The splitting can be expressed through the parameter a 
describing the effective mass polaron shift 

A 2 1fi , (h ,2 [ dq |<2|e^|l>| 2 2 m c e 4 /2(fte*) 2 

A = 167ra/ L o("-w L o) / TTTvT 2 > a = £ ■ ( 18 ) 



Here we introduced the length Zlo = -\/?i/2m c WLOj "ic is the electron effective mass. The parameter a for materials 
with a relatively weak polarity is small. For instance, a = 0.39 for CdTe with partly ionic bonding. Suppose now that 
we can achieve the resonant condition £ 2 — £1 = hco^Q changing the intcrlcvel spacing £ 2 — E\. If the states 2 and 1 arc 
the spin up and spin down states the resonant condition can be reached by adjusting external magnetic field. Since 
the level splitting is proportional to a matrix element < 2|e lqr |l > we see that the phonons can lead to spin flips only 
provided the states 2 and 1 are not the eigenfunctions of spin operators s 2 and s z . For this reason, we must include 
into the Hamiltonian the spin-orbit interaction. We consider the spin-orbit interaction in the Rashba model [19] 

H R = ^[<rp]n. (19) 

Here n is a unit vector perpendicular to the quantum well plane. This interaction is due to the structure inversion 
asymmetry. Parameter an is of the order of 1CP 9 eV-cm. 

There could be another spin-orbit interaction term that is due to the bulk inversion asymmetry. The corresponding 
3D spin-orbit Dresselhaus Hamiltonian [2(| in the crystal principal axes reads 

H D = 5((tP). (20) 
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FIG. 1: Level splitting 

Here Ti P x = p y p x p y — PzPxPz an d other components of P can be obtained by cyclic permutations. In 2D case this 
Hamiltonian takes the form (we omit the terms cubic in p) 

H D = -T-{VyPy ~ <JxPx), (21) 

where ao — 8 < p\ > /% and < p\ > is averaged over the transverse motion of the electron. The parameter ao can 
be estimated as 10 -10 eV-cm. 

Being interested only in the possibility of the line splitting in the optical reflection (transmission) experiments 
with quantum wells explicit calculations for the spin-orbit interaction of the Rashba form will be presented since in 
many semiconductor nanostructures the Rashba interaction is stronger than the Dresselhaus one. However, it can 
be shown that the Dresselhaus term in the form (|2"Tj) does not essentially differ from the Rashba term, so that to 
take into account the Dresselhaus term one should simply replace the constant an by an (this will be sufficient for 
estimations). Indeed, one can show that the Dresselhaus term can be obtained from the last term in Eq. (j50|) below 
simply replacing cir by o>d and a by — ia. 

III. DEEP QUANTUM WELL IN TRANSVERSE MAGNETIC FIELD 

Let x. y be parallel to the quantum well plane, z axis being perpendicular to the plane of the well. Further on we will 
consider the simplest case of an infinitely deep well. We assume that magnetic field B is along z axis (perpendicular 
to the plane of the well) and choose the gauge A = B(0, x, 0). In wide gap semiconductors the conduction and valence 
bands can be considered separately. In the zinc blende structures the conduction band Hamiltonian near the point 
T 6 is 

H = H +H R , (22) 

H q = ^(-iV + -^a) 2 + U + H Z . (23) 

2m c V nc J 

Here we use the basis of Ss± (where S is the S'-type Bloch amplitude and s_, s+ are the two spin functions); U is the 
confining potential of the quantum well. We write the Zeeman Hamiltonian as 

Hz = ^B<Jz9cB. (24) 

Intending to consider semimagnetic semiconductors Cdi-^Mn^Te we will incorporate into the Hamiltonian the ex- 
change Heisenberg interaction of the conduction band electrons with Mn ions 

ffce--^Jcc(r-R„)Sf 1 s, (25) 

where J cc (r — R n ) is the exchange integral of the electron with the Mn ion localized at R„ site, the sum runs over all 
the Mn ions. We will use the mean-field approximation inserting the mean value of Mn spin in z direction < S^f- n > 
instead of the corresponding operator and ascribing spin x < S^ n > to every crystal site. In this approximation the 
exchange Hamiltonian can be rewritten in the form 

H cc = -x< Sf n >N < S\J cc (r)\S >s z = -2hu c V c s z , (26) 

where No is the density of unit cells and < S\ J ce (r)|iS > is the exchange integral (that is assumed to be positive). Here 
for convenience we factor out the cyclotron frequency oj c = eB /m c c. The introduced quantity V c for the conduction 
band turns out to be negative and rather large. It can be written as 

v - * s «?Mn . N < S\J CC (*)\S > 

Vc — X <. D z > — . 

2nuj r 



The induced Mn ion spin can be written as 

<Sf" =-B s ({), C 
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k B T ' 
where Bg(£) is the Brillouin function 



For S = 5/2 



e s (0 = ^c„ tll (H£±V)-l c „ tll (i). 
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<7Mn = 2, /^s is the Bohr magneton, S = 5/2 is the spin of a manganese atom. Therefore, we see that g c in Eq. (|24|) 
must be understood as g zz — 4V C . Since N < S\J ce (r)\S >— 0.22 eV and huj c ~ lmeV we get that g c ~ 50. 

Eigenfunctions of TJo as functions of y can be chosen as plane waves e ikyV / y/Ly~. As functions of z they are the 
eigenfunctions Xi( z ) °f an infinitely deep one-dimensional well with associated eigenvalues £j. Thus one can rewrite 
-ff as 

ffo = S '-2^^ + -p( aJ -^ + ^ (28) 

Here the position of the center of oscillator xq — —k y hc/eB depends on the quasimomentum along y direction hk y 
(the motion along y axis is free) . The Rashba Hamiltonian in magnetic field is 

_ / d/dx + k y + xjll \ 

HR - aR {-d/dx + k y + x/P c o )■ (29) 



Here we have introduced the magnetic length l c = ^Jch/eB. Introducing Bose operators according to d/dx = 
(a — a t )/(v / 2Z c ), x — x = l c (a + a^)/(\/2) we get 

H = e l + Ku c {a\ + 1/2) + H z + ^!* ( ° J V (30) 

The Rashba term does not change the ground state (po(x — Xq) and its energy is Eq — Si + tku c /2 — /iBBg c /2. Other 
eigenfunctions of H are 

*<+ = ( ™ u :? n t\ x :l ) (3D 



smu n ip n+ i(x- Xq) 



with eigenvalues 



and 



with eigenvalues 



„ = -", - h,:An + I i \ ( hUJc VB 9c B \ 2 + 2 ^ (n + ^ ^ 



, -sinun^-xo) , (33) 

* cosu„<^ n+ i(a;-a;o) ' 



where 



:„ :,,/,,,(. U w('^-^^y + 2^|(n + l), (34) 



tan2u n = 2y2- 



l c [iBg c B ~ fiuj c 
Here <p n are the oscillator functions of x — xq 



<p n (x-xo) = ^tj— ===--^exp[-(x-a;o) 2 /2^] ff n [(a: - x )/l c ] (35) 

7T ' V2"n! C' 



where H n [x) are the Hermite polynomials. Thus we arrive at two groups of levels separated by (large) energy 
Hb9cB. Both groups consist of sublevels that are nearly equidistant (if one neglects the spin-orbit contribution to the 
energy) separated by the cyclotron energy Tilu c . The minimal energy in the first group is Eq, in the second group it is 
£0+ = £0 + Hb9cB- 

In what follows we restrict ourselves with the phonon induced transitions — > 0+ and n ► (n + 1) + . As we will 

see below under the realistic conditions the estimations show that the level splitting due to electron-phonon interaction 
turns out to be small as compared to the cyclotron energy, therefore it is sufficient to consider each pair of states 
separately. Indeed, for typical magnetic fields of the order of several tesla the magnetic length l c = y/hc/eB ~ 10 
nm, and the cyclotron energy Tilu c — h /mj 2 . ~ 10 -2 eV while the level splitting is of the order TlAuj ~ 10 -3 eV 
(see the estimations at the end of this section). Let us now calculate the matrix element between the ground state 
of the first group and the 0+ state of the second one 

\<i,0+,k v \e^" +i9 "\i,0,k v -g y > | 2 = sin 2 Uo | <t|e**"|*> | 2 ^expf-^Y (36) 

or taking into account that uq <C 1 

£2 _ „„ ,^LO rt , N 2„,2 j- fie 

where 



A 2 0+fi = Aa^(hco h o) 2 u 2 fl^) (37) 



f(j)= l °J dz 1 dz 2X 2 (z 1 )xU^) J dq ± ^e*p(-l 2 c ql/2-q ± \z 1 -z 2 \). (38) 

Transitions from n— to (n + 1)+ states are resonant too; for these transitions we have (omitting k y and k y — q y ) 

I < i, (n + 1) + \e iq * x+iq ' z \i, n->\ 2 = 2 ( - aR - ) I < i\e^ z \i > \ 2 

\lclJ-B9cB J 

x Me-iJi/a [L) l+1 {l 2 c q 2 j2) - L^i/2)] 2 , 

where L"(x) are the Laguerre polynomials defined as in Ref.[21]. Here we have taken into account that sinw n ~ u n ~ 
\j2(n + Vjolr/IcUbB. Since L^\{x) = L" +1 (x) — L"(x) this expression can be simplified 

2 ;2„2 



I < i, (n + 1) + |e*-'-*"*|i, n- > \ 2 = 2 ( -r-^—) \ < i\e*"\i > ^e' 1 ^ 2 [L n+1 {l 2 c q 2 j2)\ 

\IcHb9cB J 2 



We have 



. 2 
,2 _ q 'LO, , 2 / a R 



A ? „ +1)+ ,»- = ^ (R - 0) 'lcSsJ '•[!• (39) 



'l c \ f f°° l 2 q 2 

fn[j-j=lc dz l dz 2 X 2 l {z 1 )x 2 l {z2) J dq±_ -^~ 



I 2 n 2 



( 



li<f ± /2-qj_\zi-z 2 \ 



(40) 



Here we give an estimation for Aq + assuming that the transverse motion is described by the wave function Xi( z ) = 
y/2/L sin (nz/L) [see Eq.(|57|l] 

^.0^<^( ; ^)>/(t), (41, 

where f(x) = v / 2~7r/4 for x 3> 1 and f(x) = 3x/2 for a; <C 1 (see FigEJ. Taking into account the CdTe parameters, 
namely the longitudinal optical phonon frequency cjlo = 3.22 • 10 13 s _1 (246 K), the susceptibilities eo = 10.3 and 
£00 = 6.9, the effective electron mass m c — OAttiq (mo is the free electron mass) wc sec that the line splitting is 
Au> = A/U ~ an\fa/hl c ~ 10 11 s _1 . For the Dresselhaus interaction the line splitting would be ccd/cxr times 
smaller. 

IV. VALENCE BAND 

In the structures having zinc blende symmetry the valence Tg band is described by the Luttinger Hamiltonian 

H = H Q + H(k z ), (42) 
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FIG. 2: Functions f(x), fo(x) and fi(x). The second function corresponds to A^ + _ and saturates at x 3> 1 reaching the 

value 7\/27r/16 while the third one corresponds to the transition 1 > 2+ and also saturates at x 3> 1 reaching the value 

145V21F/256 



where we separate the part H(k z ) depending on k z , 



Ho = -^ {(71 + I12M + k 2 y ) - 2 l2 {J 2 x kl + J 2 k 2 y ) - A l3 {J x J y }{k x k y } + 2^JB j , (43) 



H(k z 



— | (7i + ifl2)kl - 2 l2 J 2 k 2 z - ^ 3 ({J x J z }{k x k z } + {J y J z }{k y k z }) 



(44) 



Here 71, 72, 73, k are material parameters, J is the operator of angular momentum J — 3/2, the symmetrized products 
are defined according to 



{AB} 



AB + BA 



(45) 



We add to the valence band Hamiltonian the exchange Heisenberg interaction of the valence band electrons with Mn 
ions 



ve — / J Jye{Y XV rl Jo n S 



(46) 



where J ve (v — R„) is the exchange integral of a valence band electron with a Mn ion. 
The wave function can be written as 



* = 5^F i (r)« i (r), 



(47) 



where ui (r) are the four degenerate states at the top of the valence band [22 



"±3/2 = T-7={X±iY)s±, u± 1/2 = — = 



T-^(X±iY)s T + V2Zs± 
v2 



(48) 



It is easily seen that in this basis the spin operator s z — u z /2 is also diagonal and is related to the J z operator by 
s z = J 2 /3, therefore we can rewrite the exchange Hamiltonian as 



1 



H e = ~x< 5f n > ZV < X\J ve (r)\X > -J z = -2hu, c0 V v J z 



(49) 



Here for convenience we factor out the cyclotron frequency lu c q — eB/m,QC anticipating its appearance in the following 
formulae. The introduced quantity V v for the valence band turns out to be positive and rather large. It can be 
estimated as 



V v = x < S? n > 



M _ N < X\J ve (r)\X > 



6huj r 



Here the exchange integral for the valence band < X\J ve (r)\X > is negative. 



At the typical magnetic fields of the order of several tesla the magnetic length l c = y/hc/eB ~ 10 nm, and the 
cyclotron energy hu> cQ = h 2 /m Q l 2 c ~ 10" 3 eV, while iV < X\J ve (r)\X >= -0.88 eV. Therefore V v > 1. 
We again choose the gauge A = B(0, x, 0) and introduce the operators a, a) according to 

i 1 

k x = --j=j-{a-a)), k v = -—(a + a ! ). (50) 

Replacing also the operators J x , J y by J± = J x ± iJ y we get 

Ho = -hu cQ |[ 7 i - |72 + l2J 2 z ]{a^a + 1/2) + ^(J 2 _ + J 2 + )[a 2 + (at) 2 ] + ^(J 2 - J 2 _)[a 2 - (a^) 2 ] + ^ 2 kJb} , 

(51) 

H{k z ) = -huj c0 | aJ+(J z + -) - o)J-{J z - 1/2) 1 iV2 l3 (l c k z ) - haj c0 - j 71 + 7^2 - 272 J? \ 



l2~2 l2 Ji\{l c k z y. (52) 



Further on we will use the spherical approximation, i.e. we set 72 =73- We get 



2g v J z + (71 - 772 + 72 J 2 z )(a)a + 1/2) + 1 ( 7l + ^- l2 - 2 72 J 2 \ (l c k z ) 2 (53) 



hu> c o * w 4'" " Z,K ' ' 2 {' L 2 

+yW + J-(a ] f] + {aJ+(J, + 1/2) - a)j-{J z - 1/2)} i^ l2 {l c k z ) 

where we have introduced the effective g factor in the valence band g v = hn/2 + V v . 
Due to large values of the exchange Hamiltonian g v we can omit the last two terms, i.e. 

V = {aJ+(J z + 1/2) - ^ J-(J Z - 1/2)} V2 l2 l c — + y(J+a 2 + J-(a f ) 2 )- (54) 

in Eq. (|53[) that sufficiently simplifies the problem. The reason of such a separation of the Hamiltonian is rather 
obvious, the Hamiltonian V leads to transitions changing both spin and Landau numbers and can be taken into 
account as a perturbation. In this approximation the levels can be considered independently and we have for the top 
heavy and light hole series of levels (in the hole representation) 

EUkn,n v =E 9 - 3**** + foco Ammh (n + 1/2) + — - (55) 

^-3/2 = ^n(a; - a;ofe H )x„„(^)--=u_3/2 (56) 



j-,(lh) pt , t m o(3mi + m h ) TT 2 fi 2 n 2 v 

E -l/2 nn v = E 9- huJ cOg V + ^UJ c0 (n + 1/2) + — (57) 

i/A,n,n v » Amirrih 2miL z 

Oh) (Jkyy 

W-l/2 = ^"(^ ~ a; 0fe H )Xn l ,( z )yy= M -l/2 (58) 

Here i? 9 is the gap, mi(mh) are the light (heavy) hole masses, n v is the quantization number of transverse motion and 
Xn„(z) is the corresponding wave function. We take into account that the 71,72 parameters are related to effective 
masses by 71 = m (mh + mi)/2mhmi and 72 = m (mh — mi)/imhmi. 

In this zeroth approximation phonons can not induce transitions between these states. In the next approximation 

of perturbation theory with respect to V these states are mixed and we get for the top heavy hole state 4>_ 3 / 2 

^-3/2 1 = Vo( x - x 0ky )xi(z)—f=u_ 3/2 . (59) 

V L v 

For the light hole top state we have 

,(ih) e zk y y ( 4 l2 huj cQ {l c /L) 1 

Vli/2,0,1 = ~7f= \ Vo( x ~ x ok y )Xi(z)u_ 1/2 + -— — — tpi(x - x 0ky )x2{z)u^ 3/2 \ . (60) 



L y I 3 £-1/2,0,1 — ■£'-3/2,1, 

Now it is obvious that phonon can induce transitions between these states. Suppose that sweeping the magnetic field 
we achieve the hole-phonon resonance condition between the states described by Eq. ([59]) and Eq. ([60]) 

-£■-1/2,0,1 — -£-3/2,0,1 = ftWLO 



« >r 



„ % mo(m h -mi) ( i?ft \ 

2g v huj c0 ftcJco - 2 — = ft^LO- 

4mhmi \ m L z ) 

For the value A_ x 12,-3/2 = A v describing the splitting in the valence band we get at the resonant condition 

where 

fl c \ f f 00 q 2 l 2 

^ V \~z) =lc l dzidz2 ^ Zl ^ 2 ^ Xl ^ Z2 >^ Z2 > J d Q^^ £ exp(-llq' 2 ± /2-q ± \z 1 - z 2 \). (61) 

f v (x) =x, x < 1, /„(x) = (10/97r 2 x), a; > 1. 

Let us compare the SMPR splittings in the conduction and valence bands. We evaluate 



3/4 / , / .. , , r \ 1 / 2 

Lh 2 /2m L 2 J \g c ) \m Q {m h -mi)) \f v (2y/g v /g c l c /L) 



A 0+ ,o / u R \ fg v \ ' ( -±m h rrn \ I f(h/L 



m 



and see that the splitting in the conduction band is bigger than in the valence band and is determined by the parameter 
g v / ' g c . Here l c is the magnetic length for magnetic fields required to achieve the resonance condition in the conduction 
band. 

In principle, in valence band one can also write the spin-orbital term of Rashba type [23| ] 



that in magnetic field can be rewritten as 



H vR = -HJp]n, 
n 



H vR = — — {j+a + J-a + } 



This term leads to the ratio 



3/4 / \ i/a 



A 0+: o ^ /or\ fgA I f(l c /L) 

A v ^\a'J\g c J \fl(2y/jjjf e l e /L) 



where 



fv(j) =J dz 1 dz2X 2 1 (zi)xl(z2) J^ dqq A ex V {- q 2 /2-q\z l -Z2\/lc). (62) 

Although the SMPR condition are met first for the hole states as one sweeps the magnetic field the splitting in the 
valence band turns out to be much smaller than in the conduction band. This is the consequence of the smaller 
spin-phonon coupling strength for the states strongly shifted by the Zeeman energy. 

V. RESONANT REFLECTION AND TRANSMISSION 

We consider the simplest geometry where the wave ~ e lkz is incident perpendicularly to the plane of the well. 
Neglecting in the induced current the longitudinal part (this term in the induced current has a small factor uq — 
olrJ nsgcBlc ) so that we can put V • D = et,V • E = the Maxwell equation for the wave with frequency u> can be 
written as (in this section k denotes the wave vector of light) 

—E a + k 2 E a = -^l dz'U^(z,z',ij)Ep(z r ). (63) 

Here k 2 = ui 2 eb/c 2 (we neglect the difference in the background susceptibilities of the well and barriers). We have 
taken into account that the polarization operator n a/ g (here the averaging over the distances much greater than the 
lattice parameter is implied) is 

U R {z, z', w) = I dx'dy'U R (r, r' , u). (64) 
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The Green function of operator d 2 /dz 2 + k 2 obeys the equation 

-^ + kAg(z,z') = ~5(z~z') (65) 



and is given by 



g ± (z,z')=±-e ± ^-*X (66) 



For the transmission and reflection problem one should use + (z, z') function. Then the solution of Eq. (p5|) can be 
written as 

E a = E° a e ikz - f dz'dz"g+{z, z')n%{z', z",u)E (z"). (67) 

where E^ is the amplitude of the incident wave. For z > L, where L is the width of the quantum well, we can identify 
the transmitted wave as 

El = E° a e ikz - ^Le ik * f dz'dz"e- ik *'TI%{z', z", u>)Ep(z") (68) 

and the reflected one can be identified considering z < 

El = -^e- b j dz'dz"e lkz 'u^(z', z", u)E (z"). (69) 

Assuming that H R (z, z') can be factorized as H R (z, z') = II^ 1 ' (z)U^ 2 ' (z') (such a factorization is possible since below 
we will consider transitions between two fixed states with respect to transverse motion Xn v (z),Xn a (z)) we scalarly 

multiply Eq. (|67[) by Tl a (z) and integrate over z, then we get 

T =- T ^\ dzdz'g+(z,z')n R a (z,z',Lo)+E a J dze^H^iz), (70) 



where we have introduced notation 



T= \ dz'Y^ ) {z',uj)E Si {z'). 
Jo 



Solving Eq. (|70|) for T and making use of Eq. ([68]) and Eq. (|69| we get for the amplitudes of the transmitted and 
reflected waves 



4tt/ dz'dze- ik( - z - z '^Il R p (z,z',Uj) 
2ikc 2 fi - 4?r J dz' dz e lk \ z - z '\U^(z, z', u>) 



E a - [6 a(3 + 0/ .,_ 2t . — f j_, j_„, ;fc iz-"iTTR 71 Zt , 7\ I E P< ( 71 ) 



4tt f dz'dze ik ^+ z ">IL R Jz,z',u>) 
Q " 2^fcc 2 n-47r/dz'dze ^fe l ^ - ^ 'l^« 7 (z,z^^) /3 ' V j 

In the basis e± = (e x ± ie a )/v2 in our approximation only one component of II Q ^ is nonvanishing, i.e. II ++ = 2H XX . 
Due to the symmetry relations we have Tl xx = n w = iH xy = —Hl yx . Therefore, left circularly polarized incident 
wave e_ is not reflected, while for the right polarized incident wave e + we get 

4tt f dz' dz e- ik( > z - z 'W, (z.z'w) , s 

+ 2ikc 2 h-4:7rJ dz'dze ik \ z - z '\U^ + (z,z',u)) V ' 

for the transmission coefficient E* = t + E^_e + e tkz and 

47r/ dz'dze ik ^ +z '^>n R + (z,z',u)) 



2ikc 2 Ti - Air j' dz 1 dz e ik \ z - z '\ll R ; + (z, z', u) 



(74) 



for the reflection (amplitude) coefficient E r = r + E+e + e %kz . Since the propagation direction of the wave is now 
inverted the reflected wave has the left polarization. A linearly polarized incident wave will be reflected as a circularly 
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left polarized wave. In the case where the wave length 2ir/k is bigger than the well width L (i.e. kL <C 1) the 
exponential factors can be omitted. 

Let us consider the polarization operator. We can write the formal expression for the operator 

ILVr, r', W ) = -i Y, tfU Wif lAa (r') / ^ {tanh £ -±^ (G* ( e /R + W ) - < ( e /ft + W )) Gj* ( e /7i) (75) 



AiA 



-tanh 



^.(g«( £ a)-g4( £ A))g^( £ A)}, 



where 



JAiAa(r) 



7 p?) p 

— - {$ Al (r) V $ A2 (r) - (V * Al (r) *, 2 (r) } A $ Al (r)* Aa (r). 

ZTTlg 77lgC 



(76) 



Here ( &A(r) are the eigenfunction of the Hamiltonian and Ao is the vector potential of the applied static magnetic 
field. We consider the interband transitions and assume that the valence band states are occupied while the states in 
the conduction band are empty. Keeping only the resonant contribution in Eq. (|75p we get 



A C A„ J ^ % 



(77) 



The states in the valence band we can consider as unchanged by the electron-phonon interaction (since we are interested 
only in the splitting phenomenon in the conduction band) and we get for the circularly polarized wave 



r+ 



— iT(huj — E\ — TlCdLO + £«) 



(huj — e + + e v + iO)(hco — e_ + e v + iO)/h + iT(huj — e\ — tuv^o + e v ) 
where we have introduced notation 

-lit 



r = 



hujCy/Cb 



$] [ dr'dzf XcXv (r')Jl Xc (r)- 

A C A„ J 



(78) 



(79) 



This quantity can be related to the recombination rate of the transition under consideration. Here we have taken 
into account that Tl xx = Tlyy Since we are interested in the splitting phenomenon we assume that the resonance 
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FIG. 3: Interband transitions. 
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light frequency is close to the transition from the ground state in the + group in the conduction band to the ground 
state in the valence band with J z = —1/2 (see FiglSJ) [24j. Therefore, in the following formulae we set E\ — e c \ = £q, 

£2 — £c2 — £Q+, £c2 — Eel — MSSc-D, £v — — ^-1/2. 0,1 



-£?^/„„, and A = A 0+: o- 
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The specification of the transition between the states described by Eq. (|58|) and 

e ik vV 
1p c0+ = SXnc(z) (Vo(x ~ X k y )s+ + U (pi(x - X 0ky )S-) -y= (80) 

allows us to express the recombination rate explicitly 

4tt e 2 \p cv \ 2 1 



hucy/ch 6m 2 , 2irl 2 ' 



=< S\ Px \X > . (81) 



Here we have assumed | < Xt,i|x c i > | 2 = 1 for the overlapping of the transverse quantized wave function of the 
conduction and valence bands. Now we will introduce the following dimcnsionless variables: the deviation from the 
SMPR (e c2 — e c i — ^lo)/A = <5, the optical frequency 2h(ui — luq)/A = x [ cj = (e C 2 — £v)/h being the interband 
resonance frequency], and the uncertainty in the level energy position 2hT/A = 7. Then we can write the power 
reflection coefficient R = \r + \ 2 as 

R= j 2 (x + 2S) 2 

(x + S - v/TTF) 2 (i + S + VTTs 2 ) 2 + i 2 (x + 2S) 2 ' 

If the dimensionless deviation { > 1 (i.e. the deviation from the phonon resonance condition is much bigger than 
the splitting) using Vl + S 2 ~ S we see that the single line structure is restored 

l 2 
R = „ ' „ . (83) 

x 2 + 7 2 v ; 

In the case of exact electron phonon resonance that can be achieved by sweeping the magnetic field, 6 = and we get 

2 2 

y x 

R= (x 2 -l) 2 + 7 2 x 2 ' (84 ^ 

In this case the power reflection coefficients reaches its maximal value under the optical resonant conditions. For 
linearly polarized incident wave this maximal value is 1/2. 

So far we have assumed that the energy uncertainty of the level under consideration is much smaller than the 
splitting A, otherwise the level splitting can not be resolved. Indeed, we will have for the Green function instead of 

Eq.giD 

G2(e) = e + ihT 2 -e 2 + n- (A/2) 2 /(e + i%T 1 - e 1 - hcj w + M ) (85) 

provided we take this uncertainty into account. Here we have phenomenologically introduced T 2 and Ti for the 
corresponding energy levels £2 and E\. It is seen from this expression that even for e — e\ — hcu^o — we can discard 
the second term in the denominator since A <C hTx and the level does not split. The recombination rate can be 
estimated taking into account that \p cv \ 2 /2m is of the order of the Bohr energy, huj ~ E g ~ 1.6 eV. Then it is seen 
that 7 = hT/A <C 1 as is assumed in FigfJ] 

Let us consider the case of equal Ti = T 2 widths of both levels. Then we can write for the reflection coefficient 

R. Tlfe+jff+J 

(x + 6- v /l + 5 2 +7 2 ) 2 (x + 5 + y/\ + P + 7 2 ) 2 + 4 7 2 (.t + S) 2 

where we introduce a dimensionless quantity proportional to the sum of level widths j e = AKTi/A and neglect the 
level width due to recombination processes. Fig. [5] demonstrates how the increasing of the level widths smears the 
doublet structure of the reflection line. The symmetry of this doublet structure depends on the deviation from the 
spin electron- phonon resonance (Fig.[6|). 

VI. APPLICABILITY OF PERTURBATION THEORY 

In SecHIl we considered only the simplest diagram for the self energy. Now we are going to discuss the validity of 
this approximation for the spin-phonon interaction. It is easily seen that each additional phonon line in the higher 
order diagrams can bring about additional resonant denominator, therefore we should consider the series of the most 
diverging sequence of diagrams. The situation is not unique and has been encountered earlier in the polaron problem 
in the three dimensional case and first such a consequence of diagrams has been considered by Pitaevskii [25] . 



R 
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FIG. 4: Power reflection coefficient for SMPR as a function of optical frequency x for 7 = 0.04, 0.09, 0.14. Here only the natural 
level width (due to recombination processes) is taken into account. 

R 




FIG. 5: Reflection coefficient as a function of optical frequency at 5 = 0.1 for different level widths 7 e = 0.7, 1, 1.3, 1.6. 



We consider two empty states 1 and 2 with energies £1,2 £2 = £1 + ^lo- Each state is unoccupied £1,2 > /i- 
Therefore we can write for the electron Green's function 



1/— 1 v r 7 



The phonon Green function can be written as 

D(to,r 1 ,v 2 ) =J2\C a 



-iq(ri-r 2 ) 



iq(ri-r 2 ) 



a ' q V oj — Tiuilo + i0 lj + fouled — iO 



(87) 



(88) 



where \C a \ 2 — > 27re 2 7iwLo/ ( ? 2e *- We are to evaluate the Green function for the state 2. Since we consider the empty 
electron states above the chemical potential the self energy diagrams will involve Green functions of the type 

1 



£ — ui — (£1.2 — fj,) +i0 



These functions have the pole with respect to u in the upper half-plane. Therefore we keep in the phonon Green's 
function only the part having the pole with respect to lu in the lower half-plane (otherwise the integration over u> 
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FIG. 6: Reflection coefficient for different deviations from the spin electron-phonon resonance 5 — 0, 0.2, 0.4 at 7 e = 0.5. 



vanishes), i.e. 



1 



uj — w LO + i0 
The simplest electron self energy diagram (see diagram (a) in Fig[7]) has a resonant denominator 

FJ 2 ( £ ) ~ i f ^G(e - u)D{u>) = \ — -, 

J 2n e - lu lo - (ei - MJ + *° 



(89) 



when e is in the vicinity of e 2 = £i + <^lo- Diagrams with more resonances are of two types: the first type leads 
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FIG. 7: Self energy diagrams. Resonant sections are shown by vertical lines 

to corrections to the Green function (to the line 1 in the skeleton diagram (a) in Fig|TJ)) and they can be taken 
into account regarding the Green function as renormalized, the second type leads to the corrections to the electron- 
phonon vertex. Since the corrections of the first type can be taken into account perturbatively (these diagrams do 
not involve resonant denominators) we will not consider them and concentrate on the diagrams of the second type. 
Several diagrams of the last type are presented in Fig[7l The diagrams (b) and (c) involve two and three resonant 
denominators, respectively. We can draw more complicated diagrams with two resonance denominators (similar to 
diagram (d) in Fig[7])), it is now seen that the diagrams of this type can be regarded as the diagram (b) with a block 
that does not involve resonant denominators, such a block we will call a compact block. Therefore, we can write 
the integral equation for the renormalized vertex (see Fig[8]). In the Fig [9] we show that the compact block is the 
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expansion with respect to electron-phonon coupling strength, therefore we write the integral equation keeping only 
the first term in this expansion 



+i 



r(e -u,ky- q y , e, k y , q x , q z ) =< \k y - q y I e-**-—*** 1 2k v > + (90) 

/ W" < ^ " ^^^l 2 ^ " * " «i >< 2fc y -1v- q' y \e-^ x -^ z \lk y -q' y >x (91) 

x |Cy| 2 r(e - m', k y - q' y ,e, k v ,q' x , q' z ) 

X (lu' - lu lo + i0)(e - lu - ui' - (e 2 - ft) + i0)(e - u>' - [e x - P ) + iO) ' 

Let us write this equation for the specific states 1 = {n = 0,i = 1} and 2 = {n = 0+,i = 1} (see Sec lIIII) . i.e. we 
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FIG. 8: Equation for the vertex 
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FIG. 9: Block expansion 



consider the ground states with respect to orbital motion and to spatial confinement. In order to simplify the integral 
equation we introduce the function A(q±,s,u>) by relation 



T{e-u,k y -q y ,e,k y ,q x ,q z ) = e^y-<iyn) e -<il/± < lL_^ < i\ e -^-\ l> A{q^e,u) 



(93) 



where q± = .Jq% -f q? and wave vectors are dimensionless (the factor is the magnetic length). Then, using for the 
phase factors under the integral the relation 



J{lvl' x -1^1 



i>= Yl Uaj-qDe-™^'^ 



(94) 



we get 



A(q,e,w) 



i r dw'p 3 dp _ p2/2 

+ 2j (2tt)3 6 («' 



^M0(p)A(p,e,w') 



wlo + i0)(e — lj - lj' — (e 2 - n) + iO)(s - lj' - (e± — /a) + iO) ' 



(95) 



HP) = J^ J dq z \C p .^\ 2 \< l\e^^\l> 



(96) 



Now we suppose that the function A(q, e, uj) has no poles with respect to lu in the lower half complex plane and 
consider the case when the magnetic length is much bigger than the quantum well width. The last assumption leads 
to 4>(p) = Ti/pl c and we can rewrite the integral equation for A(q, £,cl>lo) — ^(?,£) as the Fredholm equation 



A{q,e) = 1 + A / d PP 2 e - p2/2 J 2 (q P )A( P , e), 
Jo 



where parameter A includes the resonant denominator 



1 u 2 (e 2 /e*l c ) 

8 e — (ei + hcuho — fi)+i0' 



(97) 



(98) 



1G 

In reality the uncertainty of the level ih/r enter the last equation instead of iO. Let us evaluate the minimum of 
r when we remain in the framework of perturbation theory and it is then sufficient to consider only the skeleton 
diagram for the self energy. With ?Lu;lo = 0.02 eV, B = 3 T, m c = O.lmo, cxr = 10~ 9 eV-cm, a = 0.39 we get 
that the perturbation scheme is valid for the relaxation times shorter than tq = 5 • 10 -10 sec. On the other hand, to 
resolve the splitting the level uncertainty must be smaller than the level splitting A ~ 5 • 10 -4 eV. This requires the 
times bigger than 10 -12 sec. Therefore, there exists a region of relaxation times ~ 10~ n sec, where the perturbation 
theory is valid and the splitting phenomena is discernable. 

Here we wish to note that in the ordinary situation of the magnetophonon phenomena one meets the opposite case 
of big A parameter in Eq. (j97|) and the integral equation should be solved. Therefore, we suppose that the theories 
taking into account only the one phonon processes described by the skeleton self energy diagram can not be considered 
as reliable. 

VII. CONCLUSIVE REMARKS 

We have considered optical manifestation of SMPR in semimagnetic semiconductors. Due to the electron-phonon 
coupling the resonant reflection and transmission line representing the interband transitions is split into two lines. 
The distance between the lines is determined by the strength of electron-phonon coupling. 

We should, however, indicate that some points have not been taken into account in our calculation. Among them 
the most important is the natural width of the phonon levels. For the optical phonons at low temperatures it is 
determined by decay of an optical phonon into two acoustic ones. 

The natural width of the electron lines is also important. It may be determined by collisions of electrons with 
acoustic phonons and with defects of the lattice, as well as by recombination. These effects result in widening of the 
lines that has been briefly discussed. Under the conditions where these effects are strong, the lines may overlap as 
has been indicated above. 

So far we have considered a situation where the equilibrium concentration of the carriers is so low that they do not 
influence the light absorption. One can conceive, however, another case of interest where, for instance, in equilibrium 
electrons (provided by donors outside the well) fill the conduction band up to the Fermi level. In such a case transitions 
between the valence band and the states of the conduction band above the Fermi level are allowed. The oscillator 
strength for these transitions may be bigger than for those treated in the present paper. One can expect that the 
width of the electron level in the conduction band should be rather small as the electrons can emit acoustic phonons 
with the energies not bigger than the spacing between the level they occupy and the Fermi level. However, one can 
expect that the width of the level in the valence band may be much bigger. Indeed, the holes can emit phonons with 
comparatively large energies as the spacing between their level and the top of the valence band can be rather large. 

Experimental observation of SMPR can provide information about the electron-phonon interaction. Its investigation 
can also provide important information concerning various contributions into spin-orbit interaction as well as the 
strength of the exchange interaction. 
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